Abstract. We have performed a combined analysis of 
X-ray and Sunyaev-Zel'dovich data in the direction of 
the distant galaxy cluster, RX J2228+2037. Fitting a l3- 
model to the high-resolution HRI data gives Tc = 103 ± 
12 h^Q kpc and /3 = 0.54 ± 0.03. The dependency of the 
Sunyaev-Zel'dovich effect with respect to the gas tem- 
perature allows us, through the additional use of the 21 
GHz data of the cluster, to determine fcsTe = 10.4 ± 

1/2 

1.8 /lyo keV. Extrapolating the gas density profile out to 
^sq ' the virial radius {Ry — riyg = 2.9 Mpc), we derived a 
O gas mass of Mg(r < = (4.0 ± 0.2) x lO" h^^^^MQ. 
Within the hypothesis of hydrostatic equilibrium, the cor- 



o 



responding extrapolated total mass for this source is: 
b '. Mtot(r < R^) = (1.8 ± 0.4) x IO^^/i-IMq, which cor- 

^ ' responds to a gas fraction of fgas — 0.22 ± 0.06 /i^g • 
> Our results on the temperature and on the cluster mass 
OO ' classify RX J2228-f 2037 among the distant, hot and very 
' massive galaxy clusters. Our work highlights the power of 
I [ the association of galaxy cluster mapping observations in 
J> ■ X-ray and the SZ effect to derive the cluster's physical 
OO ' properties, even without X-ray spectroscopy. 
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ASTRONOMY 

AND 
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1. Introduction 

The study of galaxy clusters is one of the clues in our 
understanding of the universe. Located at the nodes of 
large-scale structures, they are the object of many differ- 
ent cosmological studies. While it is important to study 
their statistical properties, it is also important to probe 
their internal physics in the context of structure forma- 
tion in the universe. As the main baryonic component 
of galaxy clusters, the intracluster gas is a rich source 
of information. It represents between 10 and 30% of the 
cluster total mass. Largely observed in X-ray through its 
Bremsstrahlung emission, it can also be detected from sub- 
millimeter to centimeter wavelengths. In the last case, 
the signal is due to the Sunyaev-Zel'dovich (SZ) ef- 
fect ( Sunyaev fc Zel'dovich 197^ . The intracluster gas is 
highly ionized because of its high thermal temperature. 
The thermal electrons accelerate the low energy photons 
of the cosmic microwave background (CMB) by inverse 
Compton scattering. The transfer of energy from the ther- 
mal electrons to the photons produces a shift in the CMB 
spectrum towards higher energies. This characteristic dis- 
tortion of the CMB spectrum is a w ell-known signature 
of a galaxy cluster ( Birkinshaw 199£ ). Since the SZ flux 
is independent of the redshift, the study of the SZ effect 
toward distant clusters is a powerful tool to probe massive 
clusters at high redshift and to constrain the cosmological 
density parameter. 

The single dish measurements of the SZ effect allows us 
to determine the absolute intensity of the effect from cen- 
timeter to submillimeter wavelengths. When performed 
on large radio telescopes, they also are the only way to 
get high resolution SZ observations. Recent SZ studies 
at millimeter wavelengths using such observing methods 
and coupled with the use of high performance bolomet- 
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ric detectors have led to significant results. Two inde- 
pendent groups, Pointecouteau et al. ( 1999 , 2001) with 
the IRAM-30m /DIABOLO instrument and Komatsu et 
al. ( [I999i |200lD with the NRO-45m/NOBA instrument, 
have shown surprising results in the direction of the same 
target, the most luminous X-ray cluster known to date, 
RX J1347-1145. Discrepancies between the SZ and the X- 
ray (ROSAT) signal distributions have been pointed out, 
as well as evidence for an internal structuring of the in- 
tracluster gas. The existence of an internal structure is 
confirmed by the Chandra observation of RX J1347-1145 
( [Allen et al. 20011 ). 

In this paper, we present the result of a combined SZ 
and X-ray analysis in the direction of RX J2228-I-2037, 
a distant galaxy cluster (z = 0.421). It was detected in 
the ROSAT All Sky Survey (RASS). An optical identifi- 
cation of the X-ray source was given by Bade et al. ( 1998 ). 
Voges et al. ( 1999| ) identified this source as a galaxy clus- 
ter. More recently it has been included in the NORAS 
galaxy cluster survey (Bohringer et al. 200C) and in the 
ROSAT Brightest Cluster sample ( [Ebeling et al. 200d| ). 
Both catalogs are based on the RASS observations in the 
[0.1-2.4] keV energy band. The characteristics reported 
for this target in these two papers are in agreement. 
The results of the NORAS catalog are: a count rate of 
0.2 cts s~^, an estimated fiux for a temperature of 5 keV 
of Fx 

minosity of Lx 

We have observed RX J2228-h2037 through its SZ ef- 
fect using the Nobeyama Radio Observatory (NRO) 45 m 
radiotelescope at 21 GHz and combined those data with 
the X-ray image obtained from the ROSAT/HRI instru- 
ment in the [0.1-2.4] keV energy band. In the following, we 
present an extensive study of this target through the com- 
bined analysis of the two sets of data. In a first step, we 
probe the spatial structure of the cluster with the X-ray 
data (see Sec. Then, we use the direct dependency of 



4.2 X 10 ergs s ^ cm ^ and an intrinsic lu- 
20.7 X 10'*'* /i^ ergs s"*. 
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the SZ flux with respect to the gas temperature to extract 
this key parameter and furthermore we deduce the gas 
mass and the total mass of RX J2228+2037 (see Sec. |). 

Throughout this paper, we used Hq ~ 
70hro km a'^ Mpc'^, = 0.3 and 17a = 0.7. In 
this cosmology an angular scale of 1 arcmin corre- 
sponds to a physical size of 328 kpc at the redshift of 
RX J2228+2037. 



2. Observations 

2.1. The X-ray data 

Apart from the RASS data, pointed observations have 
been conducted in the direction of RX J2228-I-2037. 
A total of 30 ks exposure time has been spent 
with ROSAT/HRI on this target on the 8th and the 
12th of December 1996. The pointing coordinates were 
a20oo=22''28'"36'' and (52ooo=+20°37'12". 

W e have processed the da ta using the EXSAS soft- 
ware (Zimmermann et al.l993). In order to limit the back- 
ground events, we checked each PHA individually. We only 
selected PHA displaying a significant signal for inclusion 
in the data processing. Finally, PHA from 1 to 10 were 
selected. The events were binned into 5 arcsec pixels to 
generate a raw image of the cluster. The cluster map is 
presented in Fig |l|a. 



2.2. The SZ data 

The SZ effect observations were conducted at the 
Nobeyama Radio Observatory in December 1999 using the 
45 m radio telescope facilities at 21 GHz. The 21 GHz 
(1.43 cm) detector is a double channel equiped with low- 
noise HEMT receivers coupled to 2 GHz bandwidth back- 
ends. Basically the two channels allow detection of a circu- 
larly polarized signal. However, we have used them as in- 
dependent detectors for the same signal. This redundancy 
allows us to monitor the self consistency of the signal and 
to avoid unexpected drift and systematics in the signal 
measurement. It also increases the reliability of our mea- 
surements. The beam size is given by the diffraction limit 
of the telescope at this wavelength: Ofwhm ^ 80 arcsec, 
which corresponds to a physical angular scale of 437 kpc 
at the cluster redshift. During the whole observation, the 
system temperature was approximately 135 K. 

Observing at the 45 m telescope at 21 GHz, we were 
looking for the extended SZ signal, with the goal of map- 
ping the distribution of the SZ decrement at large an- 
gular scales. Our observations covered an area of about 
6x6 arcmin^. We directly mapped the target in the RA- 
DEC coordinates referential by using a raster/scanning 
technique: a single observation is composed by a X raster 
(the scanning is done in the right ascension direction) and 
a Y raster (the scanning is done in the declination di- 
rection). Their combination is expected to cancel out the 



scanning effects. Each single map (either a X or a Y raster) 
had a basic scanlength of 360 arcsec for a total of 9 lines. 
Two lines are separated by 40 arsec in order to achieve 
a correct spatial sampling with respect to the beam size. 
We have monitored the atmospheric emission and done 
pointing verifications in the direction of point sources ly- 
ing near the cluster location. Calibration measurements 
have been done throughout the observing run in the di- 



rection of NGC7027 (Ott et al. 1994). The main steps of 



the data processing pipeline include a baseline subtrac- 
tion (a first order polynomial is fitted to the edges of each 
line of the X and the Y rasters). The signal is cleaned 
of glitches due to cosmic high energy particles and to the 
electronic spikes in the readout systems. After the calibra- 
tion of the signal, the X and the Y rasters are combined 
using a Fourier transform filtering method to remove the 
scanning noise ( Emerson fc Graeve 198^ ). Each raster is 



reprojected on a RA-DEC grid to avoid a difference in 
the central position between observing sequences. Finally, 
the SZ map is computed from the sum of the whole set of 
rasters. 

The resulting map is presented in Fig |^b. The sen- 
sitivity level we have reached is 0.5 mJy beam~^ for an 
integration time of 34 h. 



3. Spatial analysis 

3.1. X-ray analysis 

To characterize the intrinsic emission of RX J2228+2037, 
we first cleaned the field of view of contaminating sources. 
The sources were extracted using a maximum likelihood 
method with a threshold of 8, which roughly corresponds 
to a 3.5(7 detection level. A total of 27 sources were iden- 
tified above this level and then removed. The hydrogen 
column de nsity was fixed to the g alactic value 4.68 x 



10 " cm ^ (iDickey &: Lockman 19901 ). The count rate for 



the cluster at the virial radius is 0.053 ± 0.003 cts s . 

The X-ray map seems to present an extension in the 
right ascension direction. Moreover, the shape of the clus- 
ter seems to be slightly elliptical and tilted. For this rea- 
son, we have chosen to fit different models of the intraclus- 
ter gas distribution to the data. We checked four different 
models based on a modified King profile , or so-called /3- 
model ( Cavaliere fc Fusco-Femiano 1978 ) and refered to 
them as A, B, C and D. Our models also assumed that 
the intracluster medium is isothermal. 

Models A and B are 2D elliptical models and therefore 
they have been fitted to the X-ray map. Model A is an 
elliptical /3-model. The centroid of this model has been 
fixed to the position of the X-ray maximum. Model B 
is also an elliptical /3-model, but in this second case, the 
position of the centroid is left as a free parameter. We apply 
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Fig. 1. a) ROSAT/HRI map of RX J2228.6+2037. The map has been smoothed by a Gaussian filter with 15" FWHM. 
The exposure time is 30 ks. b) SZ map of RX J2228+2037 at 21GHz. The noise level on the map is 0.5mJy/beam, 
the integration time is 34h. On both map, the overploted contours draw the significance of the X-ray detection and 
respectively correspond to 1 to 6a detection levels. The two maps are centered on the X-ray maximum position: 
acen =22''28™33.3", Seen =20°37'12.2". 



the following expression of the X-ray surface brightness in 
radial coordinates (0,6): 



Sx{OA) =^jf(0)(l + (cos2(<^-$) 
^(^-*)) (^) 



-e^ sin^(0 - 



3/3+ i 



bx 



(1) 



^^^(O) is the central X-ray surface brightness and bx 
the background value. 6c is the core radius along the ma- 
jor axis of the cluster, and e is its ellipticity, defined as 
the ratio of the major axis core radius over the minor axis 
core radius. $ is the orientation angle, the angle between 
the cluster major axis and the right ascension axis. In 
model A, Sx{0), ^c, P, e, and bx are fitted simultane- 
ously. In model B, the centroid coordinates acen and Seen 
are additional free-fit parameters. The 27 contaminating 
sources have been masked to avoid any bias in the fitting 
procedure. 

For those two models, we strictly followed the 2D fit- 
ting procedure of an elliptical /3-model described by Neu- 



mann & Bohringer (1997): we smoothed the image of the 
cluster with a Gaussian filter with a = 5 arcsec . We repro- 
duced this Gaussian filtering on the 2D elliptical /?-modcl. 
We calculated the corresponding uncertainty taking into 
account the Gauss filter. The uncertainties on the best fit 
parameters were obtained via a Monte Carlo method: a 
random Poisson noise is added to the data. This "mod- 
ified" data set is fitted the same way as the real data. 



We reproduced this procedure 100 times for different re- 
alizations of the Poisson noise. We took the value of the 
standard deviation of the 100 fitted values as the uncer- 
tainty of each parameter. 

Even though this method does not allow us to give the 
exact value of the x^i it allows us to minimize it and it 
has already shown good results for CL 0016+16, A2218 or 
3C29 5 (see Neumann & Bohringer |1997| |1999| ; Neumann 
1999|) . 



Models C and D are "Classical" ID /3-models. They 
have been fitted to the radial X-ray surface brightness pro- 
file using a maximum likelihood method. Model C is cen- 
tered on the location of the X-ray maximum, and model 
D has been fitted to the profile folded from the cluster 
centroid position (acen, Seen) determined by model B. For 
a ID /3-model, the X-ray surface brightness is: 

-3/3+i 

Sx{9) = Sx{0)\l+[-] I +bx (2) 




The radial profiles have been extracted from the X-ray 
map within a 12 arcmin circular area from the position 
chosen as the cluster center. Within this area, the X-ray 



background can be considered as a constant (Ota 2000) 



The events of the 27 contaminating sources have been 
masked, so that the resulting radial profile is free from 
any contamination (see Fig ga). The two models include 
4 free parameters: Sx{0), Oe, /? and bx- 



E. Pointecouteau et al.: SZ and X-ray combined analysis of RX J2228+2037 
Table 1. Best fit parameters for the X-ray spatial analysis 



5 



Model Sx{0) 
(10-2 



Bx 

ergs s cm~^ arcmin-^ 



) (arcsec) 



(deg 



A 2.92 ±0.49 0.35 ±0.01 17.8 ± 4.6 0.54 ± 0.06 

B 2.36 ±0.28 0.36 ±0.01 26.0 ± 5.4 0.64 ± 0.08 

C 3.08 ±0.30 0.40 ±0.01 18.9 ± 2.3 0.54 ± 0.03 

D 2.74 ±0.27 0.40 ±0.02 21.9 ± 2.2 0.56 ± 0.02 



1.35 ± 0.11 
1.26 ±0.10 



-37.2 ± 11.0 
-28.1 ± 12.0 



1.04 
1.38 



Value of the reduced x ■ See text in Sec. 3.1 for model A and B. 



Our four models have been tested against the X-ray 
data. In each case the PSF effects have been taken into 
account. The results are summarized in Table ^ The best 
fit parameters are given with their la errors bars. For 
models A and C, the location of the cluster center, which 
corresponds to the position of the maximum of the X-ray 
emission, is: acen =22''28™33.3^ Seen =20°37'12.2". In 
the case of models B and D, this position is determined 
by the best fit position for the centroid of the elliptical /3- 
model in model B: acen =22'^28™32.7", (5ce„ =20°37'11.2". 
The core radii for the different models (A, B, C, D) as well 
as the P parameters agree within the 2a limit. In case 
of the ID models, model C presents the best quality fit. 
In the following, we have used model C as the reference 
model. 



3.2. SZ analysis 

Before combining the SZ data with the X-ray data, in 
order to extract the intracluster gas temperature, we de- 
cided to check whether our SZ data were consistent with 
the results obtained from the spatial analysis of the X-ray 
data. 

If the cluster gas distribution follows a modified King 
profile, then the SZ surface brightness can be expressed 
as: 



SsziO) = SsziO) 1+ hr 



bsz 



(3) 



where Ssz{0) is the central SZ surface brightness and 
bsz the background emission level. 

Unfortunately, the quality of the SZ data is not good 
enough to allow us to test an elliptical model, so that we 
limited our analysis to a symmetrical /3-model. Using this 
distribution to model the SZ signal, we have reproduced in 
it each step of the SZ data processing (convolution with 
the beam pattern, baseline subtraction, etc) in order to 
get the most accurate model. Using a maximum likelihood 
method, we left the SZ central flux, the core radius and 
the /3 parameter as free parameters and fitted the model 
over the SZ map. The best configuration parameters are 
obtained for a reduced chi-squared value of 1.1 (96 d.o.f): 
Ssz{0) = -3.9 ± 1.5 mJy beam"!, 0^ = 20.3 ± 3.1 arcsec 



and (3 = 0.88 ±0.19. The result for 6c and (3 are compatible 
with the best fit parameter values obtained from the X-ray 
data analysis within a la and a 2a limit respectively. The 
large errors (especially in the case of (3) are due to the low 
signal to noise ratio and to the lack of resolution of our 
SZ data at 21 GHz. This made the precise determination 
of the core radius and /3 very difhcult and explained the 
strong degeneracy between those parameters. 

4. Physical properties of RX J2228-|-2037 

Knowing that the angular diameter distance of 
RX J2228+2037 is d-A = 1127 Mpc, the value of the 
angular core radius, dc — 18.9 ± 2.3 arcsec, corresponds 
to a physical core radius of Vc = 103 ± 12 h^^ kpc. From 
the best fit parameters of model C, we have extracted the 
central electron density, ne(0) — 0.014:hjQ^^ cm^'^. 

4-.1. Temperature determination 

ROSAT/HRI has no spectral resolution so that the X-ray 
data do not allow us to determine the electron temper- 
ature. However, using the high luminosity value in the 
[0.1-2.4] keV energy band reported for RX J2228+2037 
in the NORAS catalog Lx = 20.7 x 10"*" /i^ ergs 
s~^( p3ohringcr ct al. 2000| ), we can expect a quite high 
temperature for the intracluster medium. In fact, we can 
use this value of the luminosity as a lower limit for the 
bolometric luminosity of the cluster to put a lower limit 
on the gas temperature. Using the Lx — relation, we 
obtained keTe > 7.7 keV (|Arnaud fc Evrard 1999| ) or 



keTe > 6.6 keV ( |Fairlcy et al. 2000| ). 

Knowing that the SZ signal is directly proportional to 
the pressure and furthermore to the gas mass integrated 
along the line of sight, we have used the spatial informa- 
tion of our SZ map to constrain the gas temperature of 
RX J2228+2037. The central SZ surface brightness can be 
expressed as follow : 

SsziO) =yfsz{i^,Te) (4) 
where fsz{v^Te) represents the SZ effect spectral 
shape. The dependency of this spectral shape with respect 
to the gas temperature should not be neglected, especially 
in our case where we expect a high value for the tempera- 




Fig. 2. a) X- ray surface brightness profile for RX J2228+2037. The best fit f3-mode\ profile is overplotted as a solid 
line (see Sec. '3A). b) SZ radial profile at 21 GHz. The data points are displayed with their associated la error bars. 
The best fit SZ profile is shown as solid line. It corresponds to the best fit model C (see Sec. ^T] ) scaled on the SZ 
data profile with the best fit value of the gas temperature (see Sec. 4.1). 



ture. This dependency is not linear and fszi^^T^.) h as to 
be computed numerically ( Pointecouteau et al. 1998 ). y is 
the Comptonization parameter: 



y 



■ or rie (0) 



f,{l)dl 



(5) 



where fcs, ?7ie, c and ar are respectively the Boltzmann 
constant, the electron mass, the speed of light and the 
Thomson cross-section. ne(0) is the electron density at 
the center of the cluster. = riyg — 2.9 Mpc is the 
virial radius of the cluster ( Bryan fc Norman 199^ ). fg{l) 
represents the gas distribution profile along the line of 
sight. 

To determine the electron temperature from the SZ 
data, we have fixed the /3-model parameters to the best 
fit values of model C, 9c = 18.9 arcsec and j3 = 0.54 (see 
Sec. ^^ ). rie(O) has also been fixed so that the temperature 
was the only remaining free parameter in the model. We 
have tested this /3-model against the SZ map for various 
temperature values. Taking into account the temperature 
effect on the SZ spectral shape, we have derived from a 
maximum likelihood analysis a value of fesTg — 10.4 ± 
l.%h\{^ keV (x^ = 1.02). This high temperature value is 
in agreement with the lower limit value suggested by the 
[0.1-2.4] keV X-ray luminosity. The corresponding central 
Comptonization parameter is y{Q) — 2.4 x lO^''. 

From this determination and the X-ray count rate es- 
timates (see Sec. 3.1), we deduced the X-ray flux and 
luminosity of the cluster in the [0.1-2.4] energy band: 



Fx = 2.7 X 10^12 ergs g-i cni-^ and La'[0.1 - 2.4] = 
16.5 X IQ^'^hjQ ergs s^^. The X-ray bolomctric luminos- 
ity has been computed from the best fit model: L^' = 
58.3 X lO''^ ergs s^^. This allows us to check our deter- 
mination of the temperature against the Lx — relation 



expectations: ksT^ = 11.9 keV (Arnaud & Evrard 1999) 
and ksTc = 9.9 keV ([Fairley et al. 2000]). Our tempera- 



ture determination fits perfectly on the Lx — T^. relation. 
This determination confirms the fact that RX J2228-I-2037 
is indeed a hot cluster. 



4-2. Point sources 

Where the determination of the spatial parameters of 
the cluster (rc and (3) is provided by the X-ray data, 
the electron temperature is constrained by the SZ data. 
Within the hypothesis of the hydrostatic equilibrium, the 
SZ flux and the electron temperature are linearly linked, 
so that any contamination in the SZ flux measurement 
will bias the temperature estimation. Such contamina- 
tion could be inferred by the presence of point sources 
in the field of view. To prevent such biases, we have in- 
vestigated for the presence of point sources in our field of 
view. We have found three radio sources in the NVSS sur- 



vey (Condon et al. 1998), that could contaminate the SZ 



mapping of the cluster. We only know their flux measure- 
ments at 1.4 GHz, therefore we are not able to determine 
their spectra in order to extrapolate their flux at 21 GHz 
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SI 


S2 


S3 


fcsTe (keV) 




Q!2000 


22*^ 28™ 26. 18" 


22''28'"26.7r 


22''28'"32.51" 






52000 


20°37'00" 


20°39'40.2" 


20°35'32.6" 






F(1.4 GHz) (mJy) 


9.0 


19.4 


7.6 






F(21.0 GHz) (mJy) 
a = -0.5'' 
a = -1" 
free^ 


2.3 
0.6 

0.0 ±3.6 


5.0 
1.3 
0.8 ±0.4 


2.0 
0.5 

1.1 ±0.3 


9.9 ± 1.8 
10.2 ± 1.8 
11.0 ± 1.6 


2.2 
1.06 
0.93 



a is the spectral index of the power law used as a model for the radio sources emission (see Sec. 4.2) 



In this case the fluxes of the three sources are free parameters. They are fitted simultaneously with the gas temperature. 



without a few hypothesis. Assuming that the radio emis- 
sion of the sources is due to synchrotron emission, we can 
adopt a power law as a model for their spectra. To date 
many different radio sources have been observed with dif- 
ferent spectral shapes. Nevertheless, most of the spectral 
indexes of models extracted from radio source analysis fit 
into the interval: —1 < < —0.5. Using those two lim- 
its we have extrapolated the contamination sources fluxes 
from 1.4 GHz to 21 GHz using a power law spectrum: 
F{i>) oc 1/°''' . In order to estimate whether the presence 
of radio sources affects the determination of Tg, we have 
added the point source contribution to the SZ model and 
performed a new temperature estimation. In a first step, 
each point source flux was fixed to the value extrapolated 
from the 1.4 GHz measurements with the two spectral in- 
dex values of —0.5 and —1. In a second step, we considered 
their fiuxes as free parameters as well as the electron tem- 
perature. The characteristics of the radio sources and the 
results on the temperature determination are summarized 
m table I 



In each configuration, the temperature determination 
is not really affected. The three temperature values are 
consistent with the value obtained when ignoring the point 
sources. The value is bad in the case of a = —0.5. On 
the contrary, the fit for a spectral index of —1.0 is satisfy- 
ing and also agrees very well with the fluxes value we de- 
termined from the SZ data (third case) within a 2cr limit. 
This can be explained by the low fluxes expected or fitted 
at 21GHz and by the relatively high noise level on the SZ 
map. Moreover, in this case the point sources produce a 
localized effect which is washed out in the temperature 
determination over the entire map. It appears that the 
influence of the point sources present in the cluster neigh- 
borhood can be neglected in the temperature extraction. 

1/2 

Therefore the value of fc^Tg = 10.4 ±1.8 h-^Q keV can be 
taken as the isothermal temperature for RX J2228±2037. 



4-. 3. Mass estimation 

The SZ surface brightness is directly proportion to the 
pressure integrated along the line of sight. A map of 
the SZ effect is a direct projected image of the cluster 
gas mass distribution if the temperature distribution is 
isothermal. Knowing the gas temperature, the projected 
gas mass is derived from the SZ flux through a linear 
transformation : Mg{e) = C{z,Te) x FsziO). Where C 
is a constant depending on the cluster redshift and on 
the gas temperature. From our 21 GHz SZ data, we are 
able to measure the gas mass up to a radius of 0.5 Mpc: 
Mgasir < 0.5Mpc) = (2.7 ± 0.6) x lO^^hfff Mq. 

Assuming the isothermality of the cluster, we can de- 
rive the gas mass and the total mass from our best fit 



model C to the X-ray data (see Sec. 3.1) 

Mgas = fj-rup J dS J ne{l)dl 



Mtot{r) 



dlnr 



dim 



(6) 



where /imp is the average weight per massive parti- 
cle and G the gravitational constant. The measured gas 
mass from the SZ map is in good agreement with the mass 
predicted by the best fit model C: Mgas{r < 0.5Mpc) = 
2.9 X lO^hfff/^ Mq. 

Finally, we can estimate the gas mass and the total 
mass at the virial radius by extrapolating the profile cor- 
responding to the best fit parameters out to the virial ra- 
dius. Using our best fit model C which parameters are 
9c = 18.9 and f3 = 0.54 (Sec. |]|), we obtained the 
following extrapolated measurements: Mgas{r < Rv) — 
(4.0 ±0.2) X IQi-^/ifo^/^ M© and Mtot{r < Rv) = (1.8 ± 
0.4) X 10"/iTO M©. 

The errors on the mass determination take into ac- 
count the uncertainties derived from the best fit parame- 
ters. However, our results using rc and (3 for the ID and 
2D profiles do not account for a possible oblate or prolate 
form in the direction of the line-of-sight. Taking these ef- 
fects into account would certainly increase the error bars. 
However, such analysis is beyond the scope of this paper. 
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The corresponding gas fraction is fg, 

-3/2 



0.22 ± 

0.06 /i^J*'^. This result is larger than the value derived 



by Grego et al. ( 2001 ) from the SZ measurements of 18 
clusters, fg — 0.12 ±0.01, but both values agree within 
the 2(7 limit. 



5. Conclusion 

In this paper we have presented a combined analysis of 
the SZ effect and X-ray mapping observations. Through 
this analysis, we have carried out an extensive study of 
the galaxy cluster RX J2228+2037. 

This work led us to the characterization of the cluster 
physical properties, among which the gas temperature is a 
key parameter. Our estimation of this parameter, fesTg = 

1/2 

10.4 ± l.Sh^Q keV, has only been possible through the 
SZ and X-ray data combination which shows the power of 
combining this sort of data. 

This target will be soon observed by the new gener- 
ation X-ray satellites, XMM-Newton and Chandra. The 
high quality spectral data that will be provided by those 
instruments should confirm and help to improve our re- 
sults. 

Our work also emphasizes how powerful the combina- 
tion of SZ and X-ray data is. It is a preview of what we 
can expect for many sources when the full sky coverage of 
the Planck Surveyor satellite from submillimeter to cen- 
timeter wavelengths will be available. Many sources in this 
catalog are expected to have faint X-ray counterparts. The 
quality of the X-ray spectral data is expected to be poor 
in such cases so that the combination with SZ data will 
be the only path to the extraction of reliable information 
about the physics of galaxy clusters. 
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